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We study the evolution of a perfect-fluid sphere coupled to a scalar radiation field. By ensuring a 
Ricci invariant regularity as a conformally flat spacetime at the central world line we find that the 
fiuid coupled to the scalar field satisfies the equation of state pc + 3pc ~ constant at the center of the 
sphere, where the energy pc density and the pressure pc do not necessarily contain the scalar field 
contribution. The fiuid can be interpreted as anisotropic and radiant because of the scalar field, but 
it becomes perfect and non radiative at the center of the sphere. These results are being currently 
considered to build up a numerical relativistic hydrodynamic solver. 



PACS numbers: 04.25.-g, 04.25. D-, 04.40.-b 



I. INTRODUCTION 

The simulation of gravitational radiation from collaps- 
ing massive sources is one of the main motivations for 
numerical relativity P, [2| • In the characteristic formu- 
lation Isaacson, Welling, and Winicour 3 explored the 
production of gravitational waves for axially symmetric 
distributions of matter within the Winicour-Tamburino 
framework Njl. Someyears later, fundamental studies in 
ID and 3D Q, Q, were followed by a 3D code able 
to treat matter by Bishop et al. Q. A relatively sim- 
pler system is a spherically symmetric distribution of a 
perfect fluid coupled with scalar radiation Q. For other 
symmetries other than spherical, the mathematical prob- 
lem has a very similar structure (To| . In these contexts 
we can take full advantage of the characteristic approach 
to treat matter and radiation. 

We recently reported an unexpected unity in the treat- 
ment of matter by the explicit use of Bondian observers 
PH to describe matter in numerical relativity • These 
observers offer an Eulerian (non-comoving) description 
(global) with the spirit of Lagrangian observers (local 
and comoving). Consequently, wc proposed the post- 
quasistatic approximation (PQSA) jJl as a test bed 
for numerical relativity. The PQSA starts from any 
interior static solution and leads to a system of ordi- 
nary differential equations for quantities evaluated at 
the boundary surface of the fluid distribution. The nu- 
merical solution of this system allows the modeling of 
self-gravitating spheres whose static limit is the origi- 
nal "seed" solution. The approach is based on the in- 
troduction of conveniently defined effective variables and 
heuristic ansatzs, whose rationale and justification be- 
come intelligible within the context of the PQSA. In the 
quasistatic approximation, the effective variables coin- 
cide with the physical variables (pressure and density); 
the method may be regarded as an iterative method with 



each consecutive step corresponding to a stronger de- 
parture from equilibrium. It has been shown that the 
effective variables are exactly the conservative and flux 
variables in the standard ADM 3 + 1 formulation [T3 |. 

Authors commonly refer to the scalar field as a model 
of matter distribution that simplifies the treatment of the 
hydrodynamic issues. This approach has been useful to 
study non linear physics and asymptotic behaviors, es- 
pecially for central regions [ij] . Scalar field models have 
been extended to realistic situations such as gravitational 
radiation However, cases where the scalar field is cou- 
pled with radiation require a different approach due to 
both the confinement effect of matter and the dispersive 
nature of radiation. When a scalar field coupled to mat- 
ter it can be easily interpreted as an anisotropic fluid. 

We develop a numerical framework to deal with matter 
coupled to scalar radiation. We perform a detailed study 
of the central world line at r = 0, leading to a confor- 
mally flat spacetime in that region. We assume that the 
radial dependence of the geometrical and physical vari- 
ables keep the same dependence as the (quasi or post- 
quasi) static variables near the coordinate-origin. Con- 
sequently, the central equation of state (CEoS) is unique 
for all evolutions; it emerges as a conserved quantity from 
the field equations. 

In what follows we write the field equations for Bon- 
dian observers when matter is coupled to scalar radia- 
tion, which makes the fluid manifestly anisotropic. The 
asymptotic approximation to the central world line is dis- 
cussed with detail in section III, showing that the regu- 
lar spacetime is conformally flat at r = 0. Using the 
PQSA, which belongs to the aforementioned spacetime, 
we look the geometrical variables near the coordinate- 
origin. With all this background we propose an asymp- 
totic expansion which straightforwardly leads to a CEoS. 
Some remarks arc finally pointed out. 
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II. FIELD EQUATIONS 



We use the radiation metric in the spherical form [llj 



= = (1 - 2TO/r)$y2, 



2 = e^tif}^du^ + 2dudr) - r^{de^ + 



2 ' sine'^ 



\ r 



% (1) 



where (3 = f3{u, r) and V = V{u, r). In spherical symme- 
try there exists a well defined notion of quasilocal energy, 
the Misner-Sharp mass function, rh{u, r) [l^ introduced 
by means of 



(2) 



which measures the energy content in the sphere of ra- 
dius r and it reduces to the Arnowitt-Deser-Misner and 
Bondi masses in the appropriate limits. 

Now, consider an energy-momentum tensor for a per- 
fect fluid and a massless scalar field. We could fol- 
low the Tamburino-Winicour formalism in particu- 
lar as applied in regular spacetimes where the foliation 
of light cones emanates from a freely falling central ob- 
server 3,0]. But we choose a different gauge, a Bondian 
one. For an observer comoving with the fluid, with radial 
velocity u) respect to a local Minkowskian observer, the 
covariant energy-momentum tensor of matter looks like 



T, 



/9 

p 

p 

p. 



(3) 



where p and p are the energy density and the pressure 
respectively. 

The energy-momentum tensor for the massless scalar 
field minimally coupled with gravity 

T*, = V^<i>V,$ - ig^.V"<i>V„<i>, (4) 
can be characterized as a radiating and anisotropic fluid 

M 

(5) 

with w^w^i = 1, l^lf^ = 0, x'^Xm = ^1; if we identify the 
four-velocity for an observer at rest in the frame of ([T]) 
by 

u'^ = (l-2m/r)"'/'e"'''5o: (6) 
and the null and the space-like vectors by 

l, = {l-2rh/rf'^^f'5l (7) 



(8) 



The scalar energy flux, the scalar energy density, the 
scalar radial pressure, and the scalar tangential pressure, 
respectively are 



(9) 



and 



(10) 



(11) 



where the comma denotes partial differentiation respect 
to the indicated coordinate. Note that Bondian observers 
can be purely Lagrangian when we deal only with radi- 
ation. The observer with four-velocity ^ is resting at 
infinity 

Thus, the Einstein-Klein-Gordon equations can be 
written as 



P + ^ P , $ , $ 
-+ p +e 

1 — CJ^ 



uop 



47rr(r — 2m) 47rr2 



^(p + p) + P*+P* = (r-2m) ^' 
1 + UJ 



27rr2 



(12) 



(13) 



(14) 



P + Pt 



in 



-2/3 



and 



+ -Z [3/3,r(l — 2TO^r) — "T-.rr], 



2(r$),„r = r-i[re^^(r - 2m)$ , 



(15) 



(16) 



In absence of matter, p = p = 0, the field equations re- 
duce to the well known hypcrsurface equations and wave 
equation for a scalar field Q • The scalar field clearly in- 
duces anisotropy in the matter distribution. Equations 
(fT^ - p5|) are formally the same as those reported many 
years ago to study anisotropic and radiating matter [l8| . 



III. ANALYSIS OF THE SYSTEM NEAR THE 
COORDINATE-ORIGIN 

Some previous investigations consider regularization 
near r = o[il,[lil,[23l,^,[23. The conditions for the 
scalar field, as a matter model, do not necessarily apply 
to distributions of matter. Depending on gauge condi- 
tions each procedure to get regular spacetimes may be 
cumbersome and tricky, even in vacuum. Initially regu- 
lar spacetimes can eventually develop singularities [ij], 
j23t . We show a simple way to construct regular space- 
times, near the coordinate-origin, when the inner space- 
time corresponds to a spherical distribution of baryonic 
matter coupled to a massless scalar field. 
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A. Ricci invariant and conformally flat regularity 

Using the Ricci invariant 

where 

X = 2(p + p) f i-^ + i±4 P*) + 4 (p + pf ) , 



and the Weyl invariant 

^ = C^^^'^C^^.^s = 48 



(18) 



we can characterize the spacetime near the origin. The 
Ricci invariant regularity imphes that the physical vari- 
ables should be regular at r = 0. This give us from the 
field equations some restrictions on m and /?, that is, 
m(u,0) = m^r(M, 0) = m_rr(M, 0) = /^^^(wjO) = 0; there- 
fore we get '^\r=o = 0, namely, the central world line 
corresponds to a conformally flat spacetime. We will see 
below that the post-quasistatic solutions belong to this 
spacetime class. 



B. Static, slow or post— quasistatic motion 

The so called PQSA is well documented and its 
possible relevance in numerical relativity has been re- 
ported recently [l^ . For an incompressible-like model 
and for a polytropic fluid distribution, we find in both 
cases the following parity in the power expansion near 
r = for each geometrical variable: to = mj,r^ -\- m^r^ + 
• • • and /3 = /?o + /32?'^ + Pir^ + • • ■ , where the odd coef- 
ficients TO,; and the even coefficients (3i are constants or 
very specific functions of time, depending on each case. 
This radial dependence is preserved up to the PQSA. 

If near r = the matter's motion is slow, it seems then 
reasonable that the physical and geometrical variables 
preserve the radial dependence. 



C. Central equation of state 

We consider the following asymptotic expansions for 
the metric functions, which correspond near r = to 



TO = m3{u)r^ + 0{r^), 



(19) 



(20) 



and for the escalar field 

$ = $o(u) + $i(u)r + $2(w)r2 + 0{r^). (21) 



o 


np (energy density) 


np (pressure) 














1 


0Ji{l32 - 7r$i) + ■7r$i$2 


ijJi{(32 — 7r<l>i) — 27r$i<3>2 





TABLE I: Coefficients for the physical variable expansions as 
a power of r near the coordinate-origin. 



Now, using the field equations we get the physical vari- 
able expansions as a function of r as shown in Table I. 

From the evolution equation for the scalar field ((T6|) 
we get 



and 



$0 = e^/'o^i, 



$1 = 



(22) 



(23) 



From the field equation pS]) 

/32 = [7r$i$2-2wi(/32-7r$2)]e2/3o^ (24) 

and from (|12p 

2/32, 



TO3 



(25) 



where the dot indicates time differentiation. This last 
equation, together with expansions showed in Table I, is 
readily integrated to obtain 



3pc = constant. 



(26) 



Po+Po- 



where Pc = Po + Po and Pc 

This CEoS is general if the proposed expansions of the 
metric variables and the scalar field near the origin, in 
powers of the radial coordinate, is general as well. It 
is easily verifiable that lesser order expansions are not 
general enough for the used metric; a higher order does 
not give additional information. The same is not true for 
the scalar field. This type of expansion has been used in 
the construction of spacetimes that eventually develop 
critical collapse [l3|. Our main motivation to specify 
these expansions was the behavior of the metric functions 
in the case of PQSA in radiation coordinates We 
supposed that near r = all the variables (geometrical 
and physical) preserve their radial dependence, even with 
an embedded scalar field. In fact, the scalar field docs 
not affect the Ricci invariant spacetime regularity which 
is conformally flat. An interesting property is that for 
eventual maximum values of the central energy density, 
the central pressure will have to reach minimum (possibly 
negative) values in order to preserve the quantity pc+'^Pc- 
Also, according to expressions ([9]) and (fTTjl the fluid is 
perfect {pf^ = ) and non radiating (e* = 0) at r = 0; 
this result is in agreement with our expectations. The 
referential value of the constant in (pS)) corresponds to the 
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static situation (w = 0), limited physically by — 1 < w < 
1. The initial values for 7713 and (32 can be determined 
by the initial-central value of the energy-density po(0)- 
We have to specify the initial conditions for /?o and cji 
to integrate the system of two equations for (32 and 7713. 
For LJi we could fit a quadratic polynomial and update 
at any time. Our gauge leads to /3o 7^ 0. Initial values for 
<f>0 5 'i'l and (f>2 will be given by fitting the initial radial 
profile for the scalar field. For the evolution we have to 
fit only <I>2, to evolve $0 a-nd $1 using equations and 

(ESI. 



D. CEoS as the strong energy condition 

It is well known that the strong energy condition reads 



(27) 



where is the Bondian four-velocity for matter 



-1/2 



l+UJ 



1/2 



"0 



1 o~ / \ 1/2 
f — 2m/ r \ 



1-^2 



(28) 



For the present case 

R^.M^U" ^ -4^ [(p + 3p) + 2 (?7"$,a)'] < 0, (29) 
which can be written explicitly at ?' = as 

R^..W''\r=o = -47r [p, + 3pc) , (30) 



that is, Pc + '^Pc > 0, which means that the CEoS satisfies 
easily this restriction. 



IV. CONCLUDING REMARKS 



We construct a regular spacetime which is conformally 
flat at r = 0; as a realization of this class of spacetime 
we have the post-quasistatic solutions, which preserve 
with evolution the radial dependence for the inner re- 
gions. Doing a PQSA-likc treatment near the origin we 
find a CEoS as a conserved quantity for any distribu- 
tion of spherical perfect fluid matter coupled with scalar 
radiation. 

We expect to use the CEoS as a convenient test to 
check the evolution where profuse radiation could be 
present. Our disclosure of a CEoS was likely possible 
due to the characteristic formulation. Moreover, Bon- 
dian observers makes it quite easy. Its use for the central 
zone, the boundary surface, and the infinity is leading 
us to a template code which we are reporting elsewhere. 
We wonder if the Bondi-Sachs gauge [26| could lead to a 
CEoS in ADM 3+1 formulation. 
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